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Abstract 



Interferometric gravitational wave detectors could measure the frequency 
sweep of a binary inspiral [characterized by its chirp mass] to high accu- 
racy. The observed chirp mass is the intrinsic chirp mass of the binary source 
multiplied by (1 + z), where z is the redshift of the source. Assuming a non- 
zero cosmological constant, we compute the expected redshift distribution of 
observed events for an advanced LIGO detector. We find that the redshift 
distribution has a robust and sizable dependence on the cosmological con- 
stant; the data from advanced LIGO detectors could provide an independent 

measurement of the cosmological constant. 
PACS numbers: 98.80.Es, 04.80.Nn 
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I. INTRODUCTION 



A non-zero cosmological constant may help solve some of the current observational puz- 
zles, most notably, the conflict between the age of globular clusters and the apparent high 
value of the Hubble constant [which suggests a younger Universe]. [ffi] A sizable cosmolog- 
ical constant can make the Universe old in spite of a high Hubble constant, although a 
non-zero cosmological constant is ugly from the theoretical viewpoint. Whatever our aes- 
thetic preferences, the value of the cosmological constant should ultimately be determined 
by observational measurements. 

Advanced LIGO detectors can expect to observe approximately 50 neutron star binary 
inspiral events per year, from distances up to 2000 Mpc, the accuracy in the measurement 
of the signal strength can be better than 10%, and the accuracy in the measurement of the 
chirp mass [which characterizes the frequency sweep of a binary inspiral] can be better than 
0.1%. || 0. The cosmological implications of gravitational wave observations of binary 
inspiral have been discussed by several authors (§-§]]. Most recently, Finn pointed out that 
the observations of binary inspirals in an interferometric gravitational wave detector, in 
terms of the distribution of observed events with signal strength and chirp mass, can be 
quite sensitive to cosmology [/J. 

Previous discussions of the cosmological implications of gravitational wave observations 
have considered the measurements of the Hubble constant Ho [Hq = 100/ikmsec _1 Mpc _1 , 
0.5 < h < 1] and the deceleration parameter go, assuming the cosmological constant to 
be zero. Even though Markovic discussed the measurement of the cosmological constant, 
he assumed its true value to be zero. In this paper, we consider the measurement of a 
non-zero cosmological constant Q^, and find that the fl\ dependence of the observed chirp 
mass spectrum is more robust than its Hq dependence. 

Markovic proposed measuring the cosmological parameters using the observed distribu- 
tion of measured luminosity distances and (estimated) redshifts ||, while Chernoff and Finn 
suggested an alternative method without using the measured luminosity distances 0. We 
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have chosen to follow the method and notations of Refs. |J0]. 

II. THE CHIRP MASS AND THE SIGNAL-TO-NOISE RATIO 

The observed chirp mass M. = Aio(l + z), where .Mo is the intrinsic chirp mass of the 
binary source, and z is the redshift of the source. Assuming a known and constant Aio, the 
chirp mass spectrum is determined by the source redshift distribution. For simplicity, we only 
consider neutron star binaries in this paper; for typical neutron star masses, Aio = 1.19 M . 
We compute the expected redshift distribution of neutron star binary events for an 
advanced LIGO detector, assuming a non-zero cosmological constant. 

Neutron star binaries (NS/NS) may one day become the "bread and butter" sources of 
LIGO style gravitational wave detectors. NS/NS merger rate at redshift z per unit observer 
time interval per unit volume is 

^ = n (l + z)\ (1) 

where ho is the local NS/NS merger rate per unit volume, (1 + z) 2 accounts for the shrinking 
of volumes with redshift (assuming constant comoving volume density of the merger rate) 
and the time dilation. 

The "best guess" local rate density from Phinney || is 

h ~ (9.9 + 0.6 h 2 ) h x 10~ 8 Mp<r 3 yr -1 ^ 10" 7 /i Mpc~ 3 yr -1 - (2) 

i.e., 3 per year at 200Mpc for h = 0.75. Ref. |J gives h — 3h 3 x 10 -8 Mpc -3 yr _1 , which 
is consistent with Eq.@ within the uncertainty of the estimates. Since the rate of star 
formation in galaxies appears to increase rapidly as one looks back to z — 0.3, it is possible 
that we underestimate the rates for LIGOs sensitive to sources at cosmological distances 
unless we consider evolutionary effects. || In this paper, we neglect evolutionary effects. 

LIGO aims to monitor the last stage of inspiral of a NS/NS binary, during which the 
gravity waves generated sweep up in frequency, over a time of about 15 minutes, from 10 
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Hz to approximately 10 3 Hz, at which point the neutron stars collide and coalesce. The 
inspiral waveforms are determined to high accuracy by only a few parameters: the masses 
and spin angular momenta of the binary components, and the initial orbital elements (i.e., 
the elements when the waves enter the detector band). As the binary's bodies spiral closer 
and closer together, the waveform increases in amplitude and sweeps up in frequency (i.e., 
undergoing a "chirp"). The shapes of the waves, i.e., the waves' harmonic content, are 
determined by the orbital eccentricity. Gravitational radiation energy losses should lead to 
highly circular binary orbits. In the Newtonian/quadrupole approximation, for a circular 
orbit, the rate at which the frequency sweeps or "chirps", df/dt, is determined solely by the 
binary's "intrinsic chirp mass" , M = (M 1 M 2 ) 3 / 5 /(M 1 + M 2 ) 1/5 , where M 1 and M 2 are the 
two bodies' masses. The number of cycles spent near a given frequency is n = f 2 (df /dt)" 1 . 
In LIGO, the frequency sweep is monitored by matched filters. The incoming noisy signal 
is cross correlated with theoretical templates. [|1(| 



For a binary inspiral source located at redshift z, the detectors measure M. = +z), 
which is referred to as the observed chirp mass. For a given detector, the signal-to-noise ratio 



is m 

5/6 

d L \1.2M 



P = 8© ^7^7-1 C(/max), (3) 



dL is our luminosity distance to the binary inspiral source. r and C(/max) depend only 
on the detector's noise power spectrum. The characteristic distance r gives an overall 
sense of the depth to which the detector can "see". For advanced LIGO, r = 355 Mpc. |7j 
< C(/max) < 1 is a dimensionless function, its argument / max is the redshifted instantaneous 
orbital frequency when the inspiral terminates because of the coalescence (or the imminence 
of coalescence, when the orbital evolution is no longer adiabatic) of the compact components; 
( reflects the overlap of the signal power with the detector bandwidth. For source redshift 
z, ( ~ 1 for 1 + z < 10 [2.8M /(Mi + M 2 )] J7|. C — 1 is a good approximation in the context 
of this paper. is the angular orientation function, it arises from the dependence of p on 
the relative orientation of the source and the detector, < O < 4. Although can not be 



measured, its probability distribution can be approximated by 0: 

I 50(4- 0) 3 /256, ifO<0<4, 
^e(©)= (4) 
0, otherwise. 

P e (0) peaks at = 1. 

Throughout this paper, in presenting our results, we use p > p = 8, r = 355 Mpc [for 
advanced LIGO], and A4q = 1.2 M Q [for typical neutron star binaries]. 



III. MAXIMUM REDSHIFT OF OBSERVED EVENTS 



The luminosity distance cIl(z) = (1 + z) 2 d A (z), where d A (z) is the angular diameter dis- 
tance. Let us denote the fractions of the critical density as Qq (matter), Q A (the cosmological 
constant), and flk (curvature). Note that Qq + + flk = 1. Then we have 



(1 + z) d A (z) 
cH, 1 



T(z) is defined to be 



sin 

T(z) 

sinh 



T(z) 



dw 



Qk < (closed), 
Q k = (flat), 
Qk > (open). 
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In a flat Universe, 
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For n A = 0, (1+*) d A (z)/(cH 



2 [l - l/VTTz\. Forfi A = 1, (l + z)d A (z)/(cHo 1 ) = z. 
Given the detector threshold in terms of the minimum signal-to- noise ratio po, the maxi- 
mum redshift of the source that the detector can "see", z max , is given by Eq.(|3p with p = p 
and = 4, i.e., 
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where we have defined 



A = 0.4733 



^A^max) 

cHo 1 
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yp J V355Mp C y \1.2M / 
Note that z max depends on h only in the combination h A(r , p , z max increases with 

hA. 

In a flat Universe, we can find approximate analytical solution to Eq.(|8|). For Q\ > 0.5, 
z max is given by 

hA{\ + hA)- 1 ' & 



y y (1) r^j 

^max — ^max — 



[1 + z) d A (z)/(cH { 
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For close to 0, z max is given by 



"max — ^max 
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Connecting the large and small approximations with an arbitrary smoothing function, 
we have 



z = z {l) 

^max *max 



1 - e~ 5Q l 



, (2) -5fi^ 
~ Sax ° 
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The above approximate formula is accurate to better than 1.2% for Q\ < 0.5 and to better 
than 0.6% for fi A > 0.5. 

IV. DISTRIBUTIONS OF OBSERVED EVENTS IN REDSHIFT AND 

SIGNAL-TO-NOISE RATIO 



The number of binary inspiral events seen by a detector on earth per source redshift 
interval per signal-to-noise interval is 



diV(> p ) • a \A ( M 2 Cdt P ( \ 



Atttiq (cH l J 



di(z) 



n 2 



l + Z 



cHo 1 } Jn (i + zy + n A + n k (i + z) 



P p (p,z), (13) 



where dt is the differential light travel time. P p {p, z) is the probability that the detector 
"sees" a source at redshift z with signal-to- noise ratio p > p , 



p p {p,z) = p e {e{ P )) 



dp 



e 
p 



Pe (9(p)) 



(14) 



where Pe(@) is given by Eq.(f4|). 

It is straightforward to find the expected distribution of observed events in the source 
redshift z and in the signal-to-noise ratio p. The distribution in z is 



P(z,> Po) 



dN(> p )/dz _ 47m (cH x )' 



N{> po) 



N(> po) 



dA{z) 
cHo 1 



l + z 



Q {1 + zf + Q A + Q k {l + z) 2 



C e (x) 



(15) 



where Cq(x) is the probability that a given detector detects a binary inspiral at redshift z 
with signal-to-noise ratio greater than p , it decreases with z. Because of Eq.(|TJp, C®{x) = 
J™dQP {Q), hence 



C e (x) 



(l + x)(4 -x) 4 /256, if0<x<4, 



0. 
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otherwise. 



d A (z) 
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hA(r ,p ,M ) 

x is the minimum angular orientation function in terms of z, po and tq. Fig.l shows the 
distribution of observed events in the source redshift z, for Q\ = 0, 0.5, 0.9, 1, and for 
h = 0.5, 0.8, assuming A = 0.4733 and a flat Universe. P(z, > p Q ) increases with z due to 
the increase in volume with z, until Cq(x) cuts off the growth at x ~ 1. Since at small z, 
Cq(x) is most sensitive to h [the minimum angular orientation function x scales with h], 
the peak location of the redshift distribution is determined by h. Q\ determines the shape 
of the peak for given h. Note that while the /i-dependence of the peak location comes in 
through the straightforward scaling of the minimum angular orientation function x in the 
combination of h A(r , p , A4o) [which depends on detector and source properties], the shape 
of the peak depends on fl A through the angular diameter distance (^(z) in a detector and 
source independent way. 



The distribution in the signal-to-noise ratio p is 



P(P, > Po) 




d A (z)] l + z 



Pp(p,z), (18) 



N(> po) Jo 



cHo 1 ] ^/n (i + z) 3 + n A + n k (i + z) 



where P p (p, z) is given Eq . (|14D with Q(p,z) given by Eq.(Q). In a flat and static universe, 
P(p,> Po) = P(p,> Po)o — 3p|]/p 4 0. The number density of events for given signal-to- 
noise ratio p decreases sharply with increasing p. Fig. 2 shows the distribution of observed 
events in the signal-to-noise ratio p relative to the distribution in a flat and static universe, 
P(p, > po)/P{p, > po)o, for Q A = 0, 0.5, 0.9, 1, and for h = 0.5, 0.8, assuming A = 0.4733 
and a flat Universe. P(p, > po)/P(p, > po)o also depends on h through the combination of 
hA(r Q ,p ,M ). 

The total number of observed events, N(> p ), is found by integrating Eq. fll3"D over z and 
p. For the local rate density of Eq.(§), iV(> po) is more sensitive to h than to However, 



note that the h dependence comes in only through the overall factor of ho (cHq 1 ) , and in 
the combination of h A(r , p , Ai ). Fig.3 shows the total event rate per year as function 
of h in a flat Universe, assuming A = 0.4733 and ho = 10~ 7 /i Mpc~ 3 yr -1 , for = [solid 
line], 0.5 [short dashed line], 0.9 [dotted line], 1 [long dashed line]. 

V. ACCURACY IN THE MEASUREMENT OF THE COSMOLOGICAL 

CONSTANT 

To measure from the chirp mass spectrum of observed neutron star binary inspiral 
events, we need to consider the spread in the intrinsic chirp mass caused by the spread in the 
neutron star masses. This effect has been discussed by Ref. [0] in detail, and it is contained 
in our parameter A. A distribution in neutron star masses distorts the observed chirp mass 
spectrum in a somewhat symmetric manner 0, while the variations in Q A and Qo distorts 
the observed chirp mass spectrum in an asymmetric manner. We neglect the spread in the 
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intrinsic chirp mass (which is likely small) in the context of this paper, then the observed 
chirp mass spectrum is the same as P(z,p > p$). 

In a flat Universe, the accuracy of measurement of the cosmological constant is roughly 
given by the Poisson noise 1/yN [N is the total number of observed events], divided by the 
"efficiency" Q with which the change in P(z, p > po) reflects the change in the cosmological 
constant. We can estimate Q to be the ratio of the percentage change in P(z,p > po) to 
the percentage change in the cosmological constant. For HA = 0.8 x 0.4733, the accuracy in 
measuring £l A is about (6-v/iV) -1 , where N is the total number of observed events. 

Without the assumption of a flat Universe, the measurement of the cosmological constant 
is complicated by the degeneracy of different cosmological models which give rise to the same 
redshift distribution of observed events. This is because P(z,p > po) depends on both Q A 
and Qq, a given P(z,p > po) will correspond to a family of related models without the 
assumption of a flat Universe. This limits the accuracy of Q A which can be measured from 
the observed P(z,p > p ), independent of the total number of observed events. 

Q and affect P(z,p > p ) mostly through the comoving distance [see Eq.(^jj)]. The 
larger z, the more differently Qq and Q A affect P{z,p > po). For given P(z,p > po) with 
known HA, a good estimate of the correlation between the uncertainties in Qq and is 

An = ^±^-An A . (19) 

For a given cosmological model (I^a, ^o) 5 its z max is given by Eq.(Q) for given hA, the above 
equation gives a family of models (f2 A , Q' ) which give rise to approximately the same z max 
and practically indistinguishable P{z,p > po). However, when we increase tq (or increase 
^4) ; -2max increases, it becomes easier to separate models with different f2 and 
Fig. 4 shows P(z,p > p = 8) for three cosmological models 

(1) tt A = 0.4, n = 0.6; 

(2) tt A = 0,n = 0.1; 

(3) tt A = l,tt = 1.35. 

9 



These models are practically indistinguishable for r = 355 Mpc (A = 0.4733), but distin- 
guishable for r = 568 Mpc (A = 0.7573). Fig. 5 shows the corresponding distributions in 
signal-to-noise ratio p. Note that in principle, it is possible to distinguish these three models 
using the additional information from the distribution in p. Increasing A by going to a 
larger r definitely helps lifting the uncertainty in the determination of ^a and fl ; i.e., the 
measurement of ^a becomes more feasible as the detector sensitivity is increased. 

Note that the particular combination of fl\ and flo, ol = ^o(l + z rna x) 2 — ^h{z max + 2), 
is measured easily and accurately. The three degenerate curves for r = 355 Mpc in Fig. 4 
have 

(1) a = 0.22, g = -0.1; 

(2) a = 0.20, g = 0.05; 

(3) a = 0.25, g = -0.325, 

where g = Q /2 — Q\ is the deceleration parameter. Clearly, a is more accurately measured 
than g ; this is because a is similar to go [a(z max = 0) = 2g ], but contains the observed 
maximum redshift z max . At r = 568 Mpc, we have a = 0.65, 0.27, and 1.38 for the three 
models respectively. 

The larger Q\, the smaller the parameter space in (Qa, ^o) which can imitate the flat 
model with f2 A + Q = 1, hence the more accurately we can determine £7a- 

VI. CONCLUSION AND DISCUSSION 

In summary, we have calculated the expected maximum source redshift z mSuX , the source 
redshift distribution P(z,> po), the signal-to-noise ratio distribution P(p,> po), and the 
total number of events per year N(p > p ), for advanced LIGO detectors in a Universe 
with nonzero cosmological constant. z mSuX , P(z, > p ), and P(p, > po) all depend on and 
Q in a fundamental way through the angular diameter distance, and they all depend on 
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h through the combination h A(r , p Q , J^A ). N(p > p ) is very sensitive to the local binary- 
merger rate ho through no (cHq 1 ^ , the value of which is quite uncertain at this time. 

The expected redshift distribution of observed events in an advanced LIGO detector 
has a robust and sizable dependence on the cosmological constant. Although the redshift 
distribution has an apparent dependence on h which is more dominant, this dependence on 
h is superficial in the sense that it always appears in the combination of h A(tq, po, Mo) [A 
is given by Eq.@]; increasing h has exactly the same effect on the redshift distribution as 

5 /6 

increasing tq or A4 , or decreasing po, by the same amount. On the other hand, the redshift 
distribution depends on an d i n a fundamental way, this dependence is detector and 
source independent. 

If we live in a flat Universe, then the cosmological constant can be determined quite 
accurately from the expected redshift distribution of observed events with a cut on the 
signal-to-noise ratio. Assuming arbitrary geometry of the Universe, the expected redshift 
distribution of observed events with a cut on the signal-to-noise ratio may correspond to a 
family of related cosmological models with different values of Q\ and f^o [the combination 
a = Qq(1 + z m ax) 2 — ^A(z max + 2) is easily and accurately measured]; this degeneracy can be 
lifted by either accurately measuring the distribution of observed events in signal-to-noise 
ratio, or by increasing the detector sensitivity. However, it may well prove more practical 
to use other astrophysical constraints on Q\ and Q Q to resolve this degeneracy. The data 
from advanced LIGO detectors should provide an independent and robust measurement of 
the cosmological constant. 
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Figure Captions 

Fig.l The distribution of observed events in the source redshift z, for f2 A = [solid line], 
0.5 [short dashed line], 0.9 [dotted line], 1 [long dashed line], and for h = 0.5, 0.8, assuming 
A = 0.4733 and a flat Universe. 

Fig. 2 The distribution of observed events in the signal-to-noise ratio p relative to the 
distribution in a flat and static universe, for f2 A = [solid line], 0.5 [short dashed line], 0.9 
[dotted line], 1 [long dashed line], and for h = 0.5, 0.8, assuming A = 0.4733 and a flat 
Universe. 

Fig. 3 The total event rate per year as function of h in a flat Universe, assuming A = 
0.4733 and h = 10~ 7 /i Mpc~ 3 yr -1 » for fi A = [solid line], 0.5 [short dashed line], 0.9 [dotted 
line], 1 [long dashed line]. 

Fig. 4 The distribution of observed events in the source redshift z, with r = 355 Mpc, 
568 Mpc (A = 0.4733,0.7573), for three cosmological models: (1) solid line: fi A = 0.4, 
tt = 0.6; (2) dashed line: fi A = 0, fi = 0.1; (3) dotted line: fi A = 1, fl = 1.35. 

Fig. 5 The distribution of observed events in the signal-to-noise ratio p relative to the 
distribution in a flat and static universe, for the models in Fig. 4, with the same line types. 



13 




z 




1 1.2 1.4 1.6 1.8 2 



L °gioP 



Figure 3 




0.5 0.6 0.7 0.8 0.9 1 

h 




z 



